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PROPHECYdf is a Monte Carlo event generator for precise simulations of the Higgs- 
boson decay H — > ZZ/WW — > 4fermions, supporting leptonic, semileptonic, and four- 
quark final states. Both electroweak and QCD corrections are included. Treating the 
intermediate gauge bosons as resonances, the calculation covers the full Higgs-boson 
mass range above, near, and below the gauge-boson pair thresholds. In this article 
we pay particular attention to the recently implemented option of PROPHECY4f to 
generate unweighted events. 



1 Introduction 

The decay of a Standard Model Higgs boson into weak-boson pairs with a subsequent de- 
cay into four fermions, H — > ZZ/WW — > 4/, plays an important role both in the Higgs 
search at the LHC pQ and in precision Higgs physics at the planned International e + e~ 
Lincar Collider (ILC). The spin and the CP properties of the Higgs boson could be verified 
upon studying angular and invariant-mass distributions [2] of the decay fermions. In order 
to match the estimated experimental precision in predictions, a Monte Carlo generator for 
H — > ZZ/WW — * 4f including radiative corrections is needed. In the past, only the elec- 
troweak 0(a) corrections to decays into on-shell gauge bosons H — * ZZ/WW [3] and some 
leading higher-order corrections were known. However, in the threshold region the on-shell 
approximation becomes unreliable. Below the gauge-boson-pair thresholds only the leading 
order was known until recently. 

PROPHECY4f 2] is a recently constructed Monte Carlo event generator for H — > 
ZZ/WW — ► 4/ that includes electroweak and QCD corrections as well as some higher- 
order improvements. Since the process with off-shell gauge bosons is consistently considered 
without any on-shell approximations, the obtained results are valid above, near, and below 
the gauge-boson pair thresholds. In this note we briefly describe the structure of the under- 
lying calculations and illustrate the new option of PROPHECY4f to generate unweighted 
events by reproducing some of the numerical results presented in Ref. [1] . 
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2 Calculational details 

The calculation of the complete electroweak 0(a) and strong 0(a s ) corrections to the pro- 
cesses H — > 4/, which includes both the corrections to the decays H — > ZZ — > 4/ and 
H — > WW — > 4/ and their interference, is described in Ref. [I] in detail. Each ingredient of 
the calculation has been worked out twice, using independent approaches as far as possible. 

For the implementation of the finite widths of the gauge bosons we use the "complex-mass 
scheme" , which was introduced in Ref. [5] for lowest-order calculations and generalized to 
the one-loop level in Ref. [6] . In this approach the W- and Z-boson masses are consistently 
considered as complex quantities, defined as the locations of the propagator poles in the 
complex plane. The scheme fully respects all relations that follow from gauge invariance. 

The one-loop amplitudes of the virtual corrections have been generated with FeynArts, 
using the two independent versions 1 [7] and 3 [8] . They have been generated and evaluated 
both in the conventional 't Hooft-Feynman gauge and in the background-field formalism 
using the conventions of Refs. [9] and [10], respectively. One version of the algebraic part of 
the calculation is based on an in-house program implemented in Mathematica, another has 
been completed with the help of FormCalc [IT]. The one-loop tensor integrals are evaluated 
as in the calculation of the corrections to e + e~ — » 4fermions [6] [12]. They are recursively 
reduced to master integrals at the numerical level. The scalar master integrals arc evaluated 
for complex masses using the methods and results of Refs. [13] , Tensor and scalar 5-point 
functions are directly expressed in terms of 4-point integrals [14j . Tensor 4-point and 3-point 
integrals are reduced to scalar integrals with the Passarino-Veltman algorithm [15] as long 
as no small Gram determinant appears in the reduction. If small Gram determinants occur, 
the alternative reduction schemes of Ref. [TB] are applied. 

Since corrections due to the self- interaction of the Higgs boson become important for large 
Higgs masses, we have included the dominant two-loop corrections to the decay H — > VV 
proportional to G^Mj| in the large-Higgs-mass limit which were calculated in Ref. [17] . 

The soft and collinear singularities appearing in the real corrections arc treated both 
in the dipole subtraction approach |18] and in the phase-space slicing method. For the 
calculation of non-collinear-safe obscrvables we use the extension of the subtraction method 
introduced in Ref. 19J. Final-state radiation off charged leptons beyond O(a), which is 
relevant if bare lepton momenta enter the event selection, is supported for weigthed events 
only. These corrections [3] are sizeable only in regions where the lowest-order prediction is 
relatively small and can amount to 4% for muons and up to about 10% for electrons. 

3 Event generation 

PROPHECY4f employs a multi-channel Monte Carlo generator similar to RacoonWW [S] 
|2"0] and Coffer^ [Jj)][2T]. The results obtained this way have been checked using the adaptive 
integration program VEGAS [22] . In its default version PROPHECY4f generates weighted 
events, which are not positive definite. 

As a new option, the program now supports the generation of unweighted events in its 
"phase-space-slicing" branch, applying a hit-and-miss algorithm similar to the one used by 
RacoonWW. Each time an unweighted event is generated, a Fortran subroutine is called 
where information about the event is provided in the format of the Les Houches Accord [23) 
(Fortran common block HEPEUP). This subroutine can be modified by the user in order to 
read out the events. 
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In the unweighting procedure also negative events occur. Although their number is re- 
duced by using only the sum of the tree- level, the virtual, and the soft endpoint contribution, 
they cannot be avoided completely. In PROPHECY4f the remaining negative events are 
treated in the same way as the positive events, i.e. they can be read out by the user in a 
subroutine. Their contribution ranges from less than a per mille to slightly more than one 
per cent of all events, depending on the Higgs-boson mass. 

The price for generating unweighted events is an increase of CPU time by about a factor 
10 2 up to some 10 3 w.r.t. weighted-event generation, depending on the chosen 4/ final state 
and the Higgs-boson mass. The results compared below arc obtained with 5 x 10 5 unweighted 
and 5 x 10 7 weighted events. The generation of these unweighted events requires about 2 
days on a AMD Optcron 252 2.6GHz CPU. However, one should keep in mind that such 
unweighted decay events could be generated once for a chosen setup and stored in a database. 
Simulations of Higgs production at the LHC or ILC could then just randomly pick events 
for the Higgs decays from the database. 



4 Numerical results 



The input parameters and the details of the setup in our numerical evaluation are provided 
in Ref. [4], where a comprehensive survey of numerical results is presented. The results 
shown in the following are obtained without applying photon recombination, i.e. invariant 
masses and angles are derived from bare lepton momenta. 

In this brief article we focus only on the decay H — > c~c + /i~/i + and show the distributions 
in the invariant masses of the decay leptons and the angle between the Z-decay planes in 
Figs. [T] and [2j respectively. These distributions play an important role in the verification 
of the discrete quantum numbers of the Higgs boson [2J. Since the radiative corrections 
significantly distort the distributions, they have to be taken into account if these observables 
are used to set bounds on non-standard couplings. Neglecting the corrections could result in 
faking new-physics effects. A detailed discussion of the corrections to these distributions can 
be found in Ref. [4] . Here we merely emphasize the agreement between the results obtained 
with weighted and unweighted events generated with PROPHECY4E 



5 Conclusions 

The generator PROPHECY4f, which simulates the Higgs decays H -> ZZ/WW -> 4/ 
including electroweak and QCD corrections at the state of the art, is extended by an option 
for the generation of unweighted events. The consistency of the new option is illustrated in 
invariant-mass and angular distributions. 



References 



[1] S. Asai et at, Eur. Phys. J. C 32S2 (2004) 19 [hep-ph/0402~254] ; S. Abdullin et ai, 
Eur. Phys. J. C 39S2 (2005) 41. 

[2] V.D. Bar ger, K.M. Cheung A. Djouadi, B.A. Knichl and P.M. Zerwas, Phys. Rev. D 49 
(1994) 79 [hep-ph/9306270| ; S. Y. Choi, D.J. Mil ler, M.M. Miihlleitner and P.M. Zerwas, 
Phys. Lett. B 553 (2003) 61 |hep-ph/02 1007 7 ; C. P. Buszello, I. Fleck, P. Marquard 



and J. J. van der Bij, Eur. Phys. J. C 32 (2004) 209 arXiv:hep-ph/0212396|. 



LCWS/ILC2007 



'"' [io- 3 ] 



0.0025 



H 



0.002 - 
0.0015 
0.001 



0.0005 - 



i — r 

weighted 
unweighted 



e e^fi \x 

— r 



l — r 



1 



M H = 170 GeV 




95 100 



Mi 



5\%] 



60 
40 



-20 
-40 
-60 
-80 



H — ► e e + fi /i H 



1 1 1 1 1 

- M H = 170 GeV 


1 1 








H^r.-' ui "' T -.":"" 


ee, weighted 




- ee, unweighted 




jU/i, weighted 




— /i /i, unweighted 




1 1 



60 65 70 



75 80 85 
Mff [GeV] 



90 95 100 



Figure 1: Corrected distribution in the invariant mass of the pair (l.h.s.) and relative 

corrections for e~e + and pairs (r.h.s.) in the decay H — > e~e + ^~/z + , obtained with 

weighted and unweighted events. 
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Figure 2: Corrected distribution in the angle between the Z — > decay planes in the 

Higgs rest frame (l.h.s.) and relative corrections (r.h.s.) in the decay H 
obtained with weighted and unweighted events. 
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